Interaction-range effects for fermions in one dimension 



Martin Hohenadler/ Stefan Wessel,^ Maria Daghofer,^ and Fakher F. Assaad^ 

^Institutfiir Theoretische Physik und Astrophysik, Universitdt Wurzburg, 97074 Wiirzhurg, Germany 
^ Institut fiir Theoretische Festkorperphysik, JARA-FIT and JARA-HPC, 
RWTH Aachen University, 52056 Aachen, Germany 
^ Institut fiir Theoretische Festkorperphysik, IFW Dresden, 01171 Dresden, Germany 

(Dated: May 10, 2012) 

Experiments on quasi-one-dimensional systems such as quantum wires and metallic chains on 
surfaces suggest the existence of electron-electron interactions of substantial range and hence physics 
beyond the Hubbard model. We therefore investigate one-dimensional, quarter-filled chains with 
a Coulomb potential with variable screening length by quantum Monte Carlo methods and exact 
diagonalization. The Luttinger liquid interaction parameter Kp decreases with increasing interaction 
strength and range. Experimentally observed values close to 1/4 require strong interactions and/or 
large screening lengths. As predicted by bosonization, we find a metal-insulator transition at Kp = 
1/4. Upon increasing the screening length, the charge and spin correlation functions reveal the 
crossover from dominant 2A:f spin correlations to dominant 4/cf charge correlations, and a strong 
enhancement of the charge velocity. In the metallic phase, the signatures of spin-charge separation 
in the single-particle spectrum, spinon and holon bands, remain robust even for rather long-ranged 
interactions. The charge-density-wave state exhibits backfolded shadow bands. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm 



I. INTRODUCTION 

The Hubbard model has served as a framework 
to study strongly correlated electrons for almost five 
decades.^ Its relative simplicity compared to more realis- 
tic models is largely based on approximating the electron- 
electron Coulomb interaction by an onsite repulsion U 
between electrons of opposite spin. The resulting Hamil- 
tonian captures many aspects of strong correlations, in- 
cluding the Mott transition at half filling. Detailed 
knowledge about the model can be obtained by combin- 
ing the Bethe ansatz with the bosonization technique.^ 
However, experiments on quasi-one-dimensional (ID) 
systems such as quantum wires, ^ carbon nanotubes,^ or 
self-organized atom chains^ fall outside the range of va- 
lidity of the Hubbard model. This is evinced by the pos- 
sibility of an insulating, charge-ordered state at quar- 
ter filling, substantial 4/cf charge correlations, or by a 
Luttinger liquid (LL) interaction parameter smaller than 
1/2. Within a quasi-lD description, these features imply 
electron-electron interactions of finite range. 

The case of one dimension is particularly interesting 
due to the breakdown of Fermi liquid theory, the impor- 
tance of collective excitations, and the emergence of spin- 
charge separation. These phenomena can be understood 
in the framework of bosonization, which provides a 
description in terms of a few nonuniversal parameters 
valid asymptotically at long wavelengths and low ener- 
gies. In particular, knowledge of these parameters fully 
characterizes the correlation functions. 

The ID Hubbard model, describing a screened, on- 
site interaction, is a Mott insulator for any C/ > at 
half filling. Away from half filling, umklapp scattering 
is not allowed and the system remains metallic. The LL 
interaction parameter takes on values 1/2 < Kp < 1, 



leading to dominant spin density wave correlations. A 
finite interaction range permits Mott or charge-density- 
wave (CDW) transitions of the Kosterlitz-Thouless type 
at other commensurate fillings n, for example at quarter 
filling in the U-V model with onsite (/7) and nearest- 
neighbor (y) repulsion.^ In contrast to the Hubbard 
model, such transitions occur at a finite critical U deter- 
mined by the condition Kp = n^. The eflFects of extended- 
range interactions depend on the details. For example, 
the intuitive picture of long-range interactions driving 
the system to strong coupling does not always apply: for 
spinless fermions, the critical interaction for the metal- 
CDW transition is larger for the 1/r potential than for a 
nearest-neighbor repulsion;^ for spinfuU fermions, a tran- 
sition seems to be absent for the unscreened potential up 
to very strong interactions.^^ 

The 1/r Coulomb potential realized in, e.g., nanotubes 
and quantum wires, represents the extreme limit of long- 
range interactions. The logarithmic divergence of its 
Fourier transform gives rise to remarkable diflFerences, 
most notably the metallic Wigner crystal (WC) state 
with quasi- long-range Akp charge correlations, ^^'^^ and 
the existence of plasmon excitations. Strictly speaking, 
the divergence only exists for infinite systems and in the 
absence of screening. Consequently, the above phenom- 
ena are absent for any large but finite interaction range, 
and the bare Coulomb potential can be regarded as a 
special point in parameter space distinct from the LL 
liquid fixed point. The 1/r potential has been studied 
analytically^^"^^ and numerically. ^'-^^'^^'^^ 

The typical experimental situation is most likely in- 
termediate between the Hubbard limit and the bare 1/r 
potential. Within bosonization, a finite interaction range 
only leads to a renormalization of the LL parameters. ^'^^ 
However, in contrast to the Hubbard model, there ex- 
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ist no analytical methods to calculate the LL parameters 
exactly for nontrivial cases. Besides, the bosonization re- 
sults rely on a linear band dispersion, and are valid only 
at low energies and long wavelengths, a limit which is 
nontrivial to achieve both in experiment and in numerical 
simulations. On the other hand, exact numerical meth- 
ods are valid at all energies and distances and permit, 
e.g., the calculation of spectral weights of excitations. 
They provide a quantitative connection to microscopic 
model parameters, and can be used to study intermedi- 
ate interaction ranges. The ID nature of the problem 
makes numerical methods particularly powerful. 

In this work we study the effect of the electron-electron 
interaction range using exact, large-scale quantum Monte 
Carlo (QMC) simulations and exact diagonalization. The 
model chosen here makes significant simplifications over 
typical experimental situations, but we believe that our 
findings are rather general. One of the key results is 
the LL interaction parameter Kp^ which allows us to 
estimate the interaction strength and range required to 
reproduce the experimentally observed values. We also 
study the evolution of static and dynamical correlation 
functions as a function of the interaction range. Impor- 
tantly, we find that spin-charge separation in the single- 
particle spectrum is robust against increasing the inter- 
action range. Our work extends previous investigations 
of spinfull and spinless lattice models, ^'^^'^""^'^^'^^'^^ and 
continuum simulations. ^^^^ 

The paper is organized in the following way. In Sec. II 
we introduce the model and discuss related previous 
work. Section III gives details of the numerical methods. 
Our results are discussed in Sec. IV. Sec. V contains the 
conclusions. The appendix provides details about the 
application of the continuous-time (CT)QMC method. 



II. MODEL 

We consider a ID chain of length L with Hamiltonian 

L/2-1 L 

H = ^e{k)nk+ V{r)Ynini+r- (1) 

k r=0 i=l 

The kinetic term contains the usual ID tight-binding 
band structure, e{k) = —2tcosk. The electron density 
operator (summed over spin cr) at wavevector k (Wan- 
nier site i) is given by hk {rii), with hia- = cl^c-^- We 
have set the lattice constant, h and /cb equal to one, and 
take t as the unit of energy. 

The interaction matrix element V{r) is defined as 

The screening length ^ permits us to interpolate be- 
tween the Hubbard model = 0, = 2V) and long- 
range Coulomb interaction = oo, V{r) ^ 1/'^]- The 



choice (2) appears more natural than gradually adding 
more and more matrix elements for increasing distances. 
The condition r < L/2 is due to the use of periodic 
boundary conditions. V{r) as defined by Eq. (2) satisfies 
V(r) ^ as r ^ oo as well as the convexity condition 
V(r-hl)+T^(r-l) > 2V{r) for r > 1. In the classical limit 
(no hopping), this guarantees a Akp CDW ground state. 
If the second condition is not met, the competition be- 
tween 2kip and 4/cf charge order can lead to enhanced 
metallic behavior or even a CDW-metal transition with 
increasing interaction range, as observed in quarter-filled 
extended Hubbard models. ^^'^^ As we show below, our 
choice of V(r) excludes such phenomena. We have also 
compared the choice of potential (2) to an Ewald summa- 
tion for the case of Fig. 10, where the cutoff is expected 
to be most relevant, but found only minor changes in the 
form of energy shifts. 

The bosonization picture for the model (1), taking into 
account the lattice, is as follows. At half filling, any 
V > produces a Mott insulator. For commensurate 
densities n away from half filling and an interaction range 
greater than or equal to the average particle spacing 1/n, 
strong enough interactions cause a CDW transition at the 
critical point Kp = n^, beyond which umklapp scatter- 
ing becomes a relevant perturbation.^ The CDW state 
is characterized by long-range 4/cf charge order. In the 
LL phase, the dominant correlations are 2A:f spin- density 
fiuctuations for Kp > 1/3, and Ak^ charge correlations 
for Kp < 1/3. For the unscreened Coulomb potential 
with divergent Fourier transform (^ = oo), we formally 
have Kp = 0, which would suggest an insulating ground 
state, in contrast to the continuum prediction of a metal- 
lic quasi- WC made by Schulz.^^ 

The existence of a metal-insulator transition at Kp = 
has been verified numerically for the U-V model and 
the U-V1-V2 model. In contrast, for lattice fermions with 
a long-range potential (more specifically, the Pariser- 
Parr-Pople model) , numerical results^^ suggest a metallic 
ground state with the properties predicted in the absence 
of umklapp scattering. This rather surprising result, 
obtained on large but finite systems, is attributed to the 
reduction of the umklapp matrix element gs due to long- 
range interactions.^ Within bosonization, there are subtle 
but important differences between spinfull and spinless 
models (concerning umklapp scattering), and between 
odd and even filling factors (e.g., n = 1/2 and n = 1/3 
are not equivalent when considering the Luther-Emery 
point). ^ These differences seem to manifest themselves 
also in numerical studies of lattice models. For example, 
whereas spinfull fermions interacting via a 1/r poten- 
tial remain metallic even for large F,^^ a metal-insulator 
transition has been observed in the spinless case,^^ with 
the critical interaction being larger than for the extended 
Hubbard model. 

For simplicity, we consider in the following only the 
case ^ < 00, so that no divergence in the Fourier trans- 
form V{q) occurs. We further focus on quarter filling 
n = 0.5, and will see below that the model (1) is then 
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either a LL (for JTp > 1/4) or a CDW insulator (for 
K, < 1/4). 

For quarter filling, n = 0.5, most of the physics of the 
model (1) (with ^ < oc) can also be captured by simpler 
U-V or U-V1-V2 models provided the convexity condi- 
tion is satisfied. In particular, these models realize the 
non-Hubbard regime Kp < 1/2, and a metal- insulator 
transition at K = 1/4. In the metallic phase, the LL 
conjecture implies that given the same LL parameters, 
the extended Hubbard models and Eq. (1) produce iden- 
tical results, albeit with different microscopic parameters. 
However, in connection with experiments, it is crucial to 
know how strong the dependence of the LL parameters 
and hence the static and dynamical correlation functions 
on the interaction range is. We will show below that 
in order to reach the same value of Kp, the U-V model 
requires much larger (and thus rather unrealistic) inter- 
actions than a model with a larger interaction range. 



III. METHODS AND OBSERVABLES 

The majority of our results were obtained from sim- 
ulations in the stochastic series expansion (SSE) rep- 
resentation with directed loop updates. ^^'^^ The inclu- 
sion of the long-range interaction terms in Eq. (2) is 
straight forward. Due to the linear scaling of comput- 
ing time with the average expansion order, this method 
permits us to study low temperatures and long chains 
(up to L = 140 here) even in the strong-coupling regime. 
We also show results obtained with the continuous-time 
QMC method. ^^'^^ The latter is restricted to weak and 
intermediate interactions due to a less favorable scal- 
ing of computer time with temperature and system size, 
and additional numerical difficulties (see the Appendix). 
Both QMC methods are exact. 

The single-particle spectral function is of particular in- 
terest in relation to photoemission results. Since the cal- 
culation of the single-particle Green's function in SSE is 
hampered by a minus-sign problem (for periodic bound- 
aries) , we instead present results from exact diagonaliza- 
tion on clusters with L = 20. 

We consider the static charge (p) and spin (cr) structure 
factors 

S,{q) = ^e"^^ {{hrho) - {nr){no)) , 

r 

S„{q) = J2^"^''iSrS^), (3) 

r 

where Sj = ^{nj^ — nj^), and the dynamical charge and 
spin structure factors 



where pq = e*^^(n^— n)/\/Z, and \i) and \j) are eigen- 
states with energies Ei and Ej. These dynamical corre- 
lation functions can be calculated in the SSE representa- 
tion at fixed particle density and for periodic boundaries 
without a sign problem. For the analytical continuation 
we have used the maximum entropy method. 

The T = single-particle spectral function reads as 



(6) 



X S 



A-(^,.) = 5:i{^i-^-)|c,_J<^))| 



X S 



where A (A+) is related to photoemission (inverse pho- 
toemission), and IV^o^""^) denotes the groundstate for the 
sector with Nq electrons and total momentum k; the cor- 
responding energy is E^^^K In order to measure energies 
relative to the Fermi energy, we show A{k,uj — p) with 

;,= [4^e + l)_^(iVe-l)j/2. 

IV. RESULTS 

Since we used three different methods, let us state here 
that the results of Figs. 1-5, 7 and 8 were obtained us- 
ing the SSE representation. Fig. 6 with the CTQMC 
method, and Figs. 9-11 by exact diagonalization. Ex- 
cept for Fig. 2(b), results are for quarter filling n = 0.5. 

A. Luttinger liquid interaction parameter 

In the metallic regime of the model (1), the knowledge 
of the LL interaction parameter Kp together with the 
bosonization results for the correlation functions provides 
a complete description of the low-energy, long- wavelength 
physics. The crossover between the Hubbard and long- 
range cases as a function of ^, and the quantitative rela- 
tion between microscopic parameters and LL parameters, 
can be studied exactly by means of numerical methods. 
The LL parameter has previously been calculated, for ex- 
ample, for spinless fermions with a 1/r potential,^ for the 
U-V model, and for the U-V1-V2 model. 

We extract Kp from SSE QMC results for the charge 
structure factor using the relation 



Kp = lim —Sp{qi) , 



(7) 



where gi = 27r/L is the smallest, nonzero wavevector for 
a given system size, and the static structure factor is de- 
fined in Eq. (3). For each V and ^, we have performed 
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FIG. 1. Finite-size scaling of the rescaled density structure 
factor 7rSp{qi) /qi, with qi = 27r/L, for IZ/t = 3, n = 0.5 and 
^ = 0.1, 0.5, 1, 2, 3, 4, 5, 10, 20 (top to bottom). Lines are 
Hnear fits, and the extrapolated value in the thermodynamic 
limit L ^ oo defines the Luttinger liquid parameter Kp. The 
extrapolation has been carried out for all data points shown 
in Fig. 2. The system sizes were L — 32, 44, 60 forn = 0.5 and 
L = 60, 100, 140 for n = 0.1. The temperature was /3t = 2L 
for V/t = 1 and /3t = L for V/t = 3, 6, 9. 



a finite-size scaling to obtain Kp. The extrapolation is 
shown for selected values of ^ in the case F/t = 3, n = 0.5 
in Fig. 1. We find that for large enough system sizes, the 
finite-size dependence is dominated by the lowest order 
1/L, and have therefore used a linear fit for the extrap- 
olation. 

Figure 2(a) shows the dependence of Kp on V/t and 
^ at quarter filling n = 0.5. The V/t = 1 results fall 
into the Hubbard regime Kp > 1/2 for all values of ^ 
shown. For a stronger interaction V/t = 3, Kp becomes 
smaller than 1/2 for ^ ^ 2, but remains larger than 1/3, 
thereby implying dominant 2/cf correlations [see Eq. (8) 
and discussion below]. At V/t = 6, the values of Kp 
span the Hubbard, non-Hubbard and dominant 4/cf (i.e., 
Kp < 1/3) regimes. For the largest ^ = 20, the LL 
parameter takes on almost exactly the critical value Kp = 
1/4 of the LL-CDW transition. The numerical results 
therefore suggest that the experimentally observed values 
of Kp ^0.25 require surprisingly large values of V/t and 
^. Finally, for V/t = 9, the system undergoes the metal- 
insulator transition for ^ ^ 3.5. Independent of V, we 
expect Kp ^ for ^ ^ oo in the thermodynamic limit, 
corresponding to the quasi- WC. A theoretical prediction, 
Kp ~ ln~^/^ ^, was made by Schulz,^^ and the numerical 
results for the charge structure factor by Fano et al}^ 
are consistent with Kp = 0. 

Figure 2(a) reveals that Kp decreases with increasing 
^, thereby bringing the system closer to the insulating 
phase. In previous work on extended Hubbard models, it 
was found that adding interactions at distances beyond 
the interparticle spacing 1/n can increase Kp and hence 
enhance the metallic character of the system. ^'^^ Simi- 
larly, in the U-V1-V2 model with U fixed, varying the 
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FIG. 2. (Color online) Luttinger liquid parameter Kp as a 
function of screening length ^. Points represent values ob- 
tained from a finite-size scaling (see Fig. 1). Lines are guides 
to the eye. (a) Results at quarter filling n = 0.5. (b) Fixed 
V/t = 3 and different fillings n. The phases are a Luttinger 
liquid (LL) for Kp > and a charge-density-wave insulator 



(CDW) for Kp < 



relative strength of Vi and V2 leads to a competition be- 
tween 2A:f and 4/cf charge fluctuations. ^^'^^ As a result, 
Kp takes on a maximum for V2 = Vi/2, where the metal- 
lic state is most stable, and it is not clear if the U-Vi- 
V2 becomes insulating at finite values of Vi and V2-^^'^^ 
The condition V{2) = V{l)/2 is also realized for the un- 
screened Coulomb potential, and numerical results sug- 
gest that the system remains metallic up to very strong 
interactions even in the presence of a lattice. ^^'^^ The 
experimentally motivated form (2), fulfilling the mono- 
tonicity and the convexity condition, favors a 4A:f CDW 
state in the limit V/t 00.^^ Similar to previous results 
for spinless fermions with a 1/r potential,^ Kp in Fig. 2 
decreases with increasing V/t. 

A common feature of the curves in Fig. 2(a) is a pro- 
nounced decrease at small values of ^, followed by a much 
slower decrease for larger ^. The numerical results indi- 
cate that the change in behavior occurs when the interac- 
tion range ^ equals the average particle spacing 1/n = 2. 
To verify this hypothesis, we compare in Fig. 2(b) the ^ 
dependence of Kp for two different densities n = 0.5 and 
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n = 0.1 ditV/t = 3. The curve forn = 0.1 indeed exhibits 
a significant ^ dependence up to much larger ^. The re- 
sults for n = 0.1 further reveal that for a given V/t, a 
smaller density requires a significantly larger interaction 
strength and/or range to reach the critical Kp = for 
the metal-insulator transition, see also Ref. 21. 

B. Charge and spin correlation functions 

For a model with SU (2) spin symmetry such as Hamil- 
tonian (1), bosonization predicts the decay of charge and 
spin correlation functions to be determined solely by the 
parameter Kp (since K^t = 1),^^ 

^''"'"'^ " + cos(2fcFx) + ^ cosi4kFx) , 

= -jdxy + ^ '^^'^^^^^^ ■ 

The 1/x^ dependence of the leading term in both chan- 
nels is familiar from Fermi liquid theory. The 2A:f 
and charge correlations decay to leading order as 
x~^~-^p and x~^^p, respectively. For the Hubbard 
model, Kp > 1/2 and 2/cf correlations dominate; tak- 
ing into account logarithmic corrections not included in 
Eq. (8), the dominant correlations in this regime are 2/cf 
spin correlations.^ For models with a nonzero interaction 
range, the 4/cf density oscillations can become dominant 
for Kp < 1/3. For the Hubbard model, subdominant 
4/cf oscillations have been observed in systems with open 
boundary conditions. 

In the opposite limit of a 1/r Coulomb potential = 
oo) with divergent Fourier transform V{q) ^ ln(l/g), 
Schulz^^ obtained 

(n^no) = ^e-"2^^cos(2i^F^) 

X 

(S^S^) = cos(2/cfx) . (9) 

Apart from the absence of the 1/x^ Fermi liquid contri- 
bution, the most notable difference is that charge cor- 
relations are dominated by an unusually slow decay of 
the 4/cf component (slower than any power law). These 
quasi-long-range 4/cf charge oscillations led to the no- 
tion of a fluctuating WC, where the wavelength A = 
27r/4A:F = l/n is the average distance between fermions. 
In contrast, the spin sector retains a power-law decay. 
These continuum results are consistent with numerical 
work.10'27 

As emphasized before, the WC results (9) rely on the 
divergence of the Fourier transform of the potential V{r). 
Such a divergence only occurs in the thermodynamic 
limit, and for ^ = oo. If either of these conditions is 
not met, the LL forms (8) can be recovered in the long- 
wavelength limit. Here we only consider large but fi- 
nite values of ^, for which a metal-insulator transition 
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FIG. 3. (Color online) Charge [(a), (b)] and spin [(c), (d)] 
structure factors Sp/a{q) for different values of the screening 
length ^, and I3t — L — 44. The key in (a) applies to (a)- 
(d). (e) Rescaled charge structure factor for y/t = 3 (/3t = 
L — 60), revealing the asymptotic approach to the small-^ 
behavior of the WC, Sp[q) ^ ^/|ln^|^/^. Here and in all 
subsequent figures, results are for quarter filling n = 0.5. 



occurs at Kp = = 1/4. The CDW state exhibits long- 
range 4/cf charge order. The closest analog of the metal- 
lic quasi- WC state in our case is therefore the metallic 
regime 1/3 > Kp > 1/4 with dominant (power-law) 4/cf 
correlations. As shown in Fig. 2, Kp < 1/3 is realized for 
V/t = 6 and large ^, and we explore the similarities to 
the WC below. 

Figure 3 shows the charge and spin structure factors 
as defined in Eq. (3). At V/t = 3 and with increasing 
^, we see a slight increase of the 4A:f = tt charge correla- 
tions, see Fig. 3(a). This effect becomes more noticeable 
for a stronger repulsion V/t = 6, as shown in Fig. 3(b). 
The inherent length scale l/n again appears in Fig. 3, 
with the results saturating on the scale of the plots for 
^ > 2. The spin structure factor [Figs. 3(c) and (d)] 
reveals an enhancement of 2/cf = 7t/2 antiferromagnetic 
correlations with increasing ^, which according to Eq. (8) 
can be related to the reduction of Kp. This enhancement 
is again more pronounced for V/t = 6 than for V/t = 3. 

Let us now turn to the long- wavelength behavior. 
For a LL we have Sp{q) ~ qKp^ whereas for the WC 
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FIG. 4. (Color online) Density-density correlations in real 
space (symbols). Here Vji — 5, /3t — L — 84, n — 0.5 and 
(a) ^ = 10, (b) ^ = 20. Lines are fits to the LL result for 
(rice no) with fitting parameters Ai [see Eq. (9)] and Kp 
determined from a (linear) finite-size extrapolation based on 
L = 44, 84. The fitting interval was (a) [15 : 35], (b) [25 : 45]. 



Sp{(i) - ^|ln(7|-i/2 (see Ref. 2). Following Ref. 10, we 
plot in Fig. 3(e) 5'p(^)| lng|-^/^/g. This quantity shows a 
logarithmic divergence at g = as long as Spicj) ~ q and 
tends to a constant as g ^ for ^ = oo.^^ Our numeri- 
cal results show that a divergence occurs throughout the 
metallic phase, and that the approach to the WC result 
is rather slow. In particular, given the finite values of 
^, the LL nature of the system reemerges eventually in 
the limit g ^ 0, although the system sizes required to 
see this effect become larger and larger. A nonlinear (at 
long wavelengths) density structure factor corresponding 
to Kp = has been observed for the 1/r potential. In 
contrast, for finite ^, Fig. 3 shows that the linear behavior 
of Sp[(i) is preserved. The long- wavelength spin structure 
factor is not affected by the interactions [Fig. 3(c) and 
(d)]; the slope in the limit q ^ ^ remains fixed, as re- 
quired by = 1 [cf. Eq. (8)]. 

Schulz^^ suggested that for a finite ^, one should be 
able to observe WC-like correlations at distances x < 
^ and LL-like correlations at x > ^. Although the 
bosonization results are only valid for large distances, 
this prediction can in principle be tested numerically. 
Figure 4 shows the density-density correlation function 
in real space. We have chosen V/t = 5, and ^ = 10 or 
^ = 20. This choice was made for the following reasons. 
First, deviations from the LL form given by Eq. (8) are 
most visible in the regime where 4A:f oscillations domi- 
nate, that is for Kp < 1/3. However, for the bosonization 
results to apply, it is important to avoid the insulating 
state expected for Kp < 1/4. Close to Kp = 1/4, pre- 
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FIG. 5. (Color online) Finite-size scaling of the amplitude of 
4kF charge correlations at fixed ^ = 10. The results reveal 
the absence of long-range order for V/t = 3, and long-range 
4A;f charge correlations for V/t = 9. The case IZ/t = 6 is on 
the metallic side but very close to the critical point, and we 
find a small but finite extrapolated value. The lines are linear 
fits. Here /3t = L and n = 0.5. 
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FIG. 6. (Color online) (a) Dynamical charge and (b) spin 
structure factor for V/t = 1 and ^ = 10. Results were ob- 
tained with the projective CTQMC method*^^ using L = 28, 
= 15 and n = 0.5. Dashed lines indicate the velocity of 
long- wavelength charge and spin excitations. 



vious work on the extended (U-V) Hubbard model has 
shown the importance of logarithmic corrections.^^ For 
the parameters chosen, we have Kp ^ 0.29 for ^ = 10 
and i^p ?^ 0.28 for ^ = 20. The results in Fig. 4 show 
dominant 4A:f correlations but no long-range order, as 
expected in the LL regime. 

Based on the idea that the LL form for {nxfio) should 
hold at distances larger than ^, we fit the numerical 
data to Eq. (8) using two fitting parameters (the 2/cf 
and 4/cf amplitudes) as well as the above values of Kp. 
The fitting intervals are chosen as [^ + 5, 35] and we used 
Pt = L = 84. Figure 4(a) shows that we indeed have good 
agreement between the fit and the QMC data at large dis- 
tances. However, for r < ^ = 10, significant deviations 
become visible. To discriminate between short-distance 
effects coming from the continuum approximation un- 
derlying Eq. (8) and genuine deviations from LL theory 
we consider ^ = 20 in Fig. 4(b). Again there is reason- 
able agreement at large distances, but clear differences at 
r < ^ = 20. Hence, keeping in mind the difficulties men- 
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FIG. 7. (Color online) Dynamical charge [(a)-(c)] and spin 
[(d)-(f)] structure factor from simulations in the SSE repre- 
sentation for y/t = 6 and different screening lengths ^. Re- 
sults are for n = 0.5, L = 44 and — 2L. The dashed lines 
indicate the velocity at long wavelengths. 



tioned above, our results are consistent with the picture 
proposed by Schulz.^^ 

As can be seen from Fig. 2, the insulating CDW phase 
can be reached for i > 3.5 and V/t = 9. The CDW state 
is characterized by long-range 4A:f charge order at T = 0, 
as formally reflected by Eq. (8) for Kp = 0, and may be 
regarded as a WC pinned to the lattice. Figure 5 shows 
the amplitude of 4A:f charge correlations divided by sys- 
tem size, i.e., Sp{AkY)/L. At fixed ^ = 10, we find that 
this quantity extrapolates to zero in the thermodynamic 
limit in the LL phase [Kp = 0.383(1), V/t = 3], and to a 
finite value in the CDW state (V/t = 9). Near the phase 
boundary, the Kosterlitz-Thouless nature of the transi- 
tion makes numerical studies difficult and we see that, 
assuming a linear scaling, Sp{AkY)/L extrapolates to a 
finite but very small value despite Kp = 0.270(3) > 1/4. 
For the unscreened Coulomb potential, Sp{AkY)/L in- 
creases logarithmically with system size, and there is no 
long-range order. 
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FIG. 8. (Color online) (a) Dynamical charge and (b) spin 
structure factor for V/t — 9 and ^ = 10, corresponding to 
the insulating CDW phase. Results were obtained in the SSE 
representation using n = 0.5, L = 44 and I3t — 2L. 



C. Dynamical charge and spin correlations 

We now discuss the dynamical spin and charge correla- 
tion functions, defined in Eq. (4), as obtained from QMC 
simulations. We begin with a rather weak interaction 
V/t = 1 and a large screening length ^ = 10. CTQMC re- 
sults for these parameters which, according to Fig. 2, fall 
into the Hubbard regime, are presented in Fig. 6. Despite 
the long-range interaction, the spectra closely resemble 
previous results for the Hubbard model, see, e.g., Ref. 40. 
In particular, the particle-hole continuum is clearly visi- 
ble in both the charge and the spin channels. As a result 
of interactions, the velocities of long- wavelength charge 
and spin excitations differ by about a factor of 2. 

To investigate larger values of V/t, we use the SSE 
representation. The latter can also be used for the pa- 
rameters of Fig. 5, but we chose the CTQMC method to 
demonstrate its applicability to models with long-range 
interactions. Taking V/t = 6, we can explore the whole 
metallic regime of the model (1) by varying the screening 
length ^. Results are shown in Fig. 7. 

We first discuss the charge sector. For ^ = 0.1, corre- 
sponding to the strong-coupling regime of the Hubbard 
model [U = 2V(0) = 12t], the results in Fig. 7(a) look 
qualitatively similar to Fig. 6(a). However, the distribu- 
tion of spectral weight over the particle-hole continuum 
is much more inhomogeneous, with pronounced excita- 
tion features along the edges. The charge velocity Vp is 
only slightly smaller than in Fig. 6(a). Upon increasing 
^, we observe a substantial increase of Vp^ as indicated 
by the dashed lines; between ^ = 0.1 and ^ = 1, 'Up in- 
creases from 1.97(2)t to 2.64(2)t. A small charge gap 
of order O.lt, which extrapolates to zero for L ^ oo in 
the LL phase, is visible in Fig. 7(c), but we can estimate 
the velocity as Vp > 3.5t. The increase of Vp reflects the 
fact that the extended interaction promotes 4A:f charge 
order, and thereby increases the stiffness of the charges 
with respect to long-wavelength excitations. This gap is 
a finite-size effect caused by the close proximity of the 
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FIG. 9. (Color online) Single-particle spectral function 
A(k^L>j — \i) for n — 0.5, Vfi — 1 and ^ = 10 from exact 
diagonalization with L — 20. We used an artificial broaden- 
ing of 0.05t. The curves marked spinon, holon and shadow 
are explained in the text. 

CDW transition. The onset of 4/cf fluctuations is also 
reflected in an incomplete but well visible softening of 
the excitations at g = 4/cf. We will see below that this 
feature develops into a Bragg peak in the CDW state. A 
plasmon excitation, one of the hallmark features of the 
1/r Coulomb potential, is not expected for finite values 
of ^, and would in general be very difficult to distinguish 
from a linear mode in numerical simulations. 

In contrast to the charge sector, the effect of ^ on the 
spin dynamics is very small. In accordance with LL the- 
ory, the velocity Vcr of long- wavelength spin excitations 
remains virtually unchanged upon increasing ^ from 0.1 
to 10 [Fig. 7(d) and (f)]. However, Vcr is strongly renor- 
malized in going from V/t = 1 [Fig. 6(b)] to V/t = 6 
[Fig. 7(d)]. At fixed V/t, the screening length hence pro- 
vides a natural way of changing the ratio of charge and 
spin energy scales, and opens a route to explore the spin- 
incoherent LL.^^ 

Figure 8 shows results for the charge and spin dynam- 
ics in the CDW phase, for V/t = 9 and ^ = 10. As 
demonstrated in Fig. 5, for these parameters, the system 
is in a CDW state with long-range 4/cf order. In addi- 
tion to a charge gap at g = 0, the charge structure factor 
has become almost perfectly symmetric with respect to 
q = 7t/2. This doubling of the unit cell results from the 
softening ai q = 4/cf, and is a typical signature of the 
CDW state. Except for a smaller velocity v^^ the spin 
structure factor in Fig. 8 is similar to the metallic regime 
(i.e., gapless), see for example Fig. 7(c). 

D. Single-particle spectral function 

The single-particle spectrum is of particular interest 
in the search for experimental realizations of LLs be- 
cause it can reveal the signatures of spin-charge sepa- 
ration (spinon and holon bands). ^^'^^ Although LL the- 



ory is a low-energy description, spin-charge separation 
may be observed up to rather high energies. For ex- 
ample, spinon and holon bands are visible over an en- 
ergy range of the order of the bandwidth in the Hubbard 
model,'^^''^'^''^^ and also experimentally for TTF-TCNQ"^^ 
and ID cuprates.^^'^^ In contrast, such clear features of 
spin-charge separation seem to be absent in recent mea- 
surements on self-organized gold chains, although the 
density of states reveals the scaling expected for a LL.^'^^ 

To understand the role of the interaction range and 
small values of Kp, we calculate the single-particle spec- 
tral function A{k,uj — ji) [Eq. (6)] for different values of 
V and ^. To simplify the interpretation of the complex 
structures, we use exact diagonalization on chains with 
L = 20 sites, and use a different graphical representation. 

Figure 9 shows the single-particle spectrum in the Hub- 
bard regime for V/t = 1 and ^ = 10. To highlight 
the spinon, holon and shadow bands previously observed 
for the Hubbard model away from half filling, ^^'^^'^^ 
we include the holon and shadow band dispersions for 
the = oo Hubbard model, — 2tcos(|A:| + k^) and 
— 2tcos(|/c| — /cf), as well as a linear spinon branch 
Vcr{k — kF) with Vcr determined from Saiq^oj). These ana- 
lytical results have well-defined corresponding excitations 
in the numerical spectra, and establish the signatures 
of spin-charge separation in the Hubbard regime of the 
phase diagram. The spectral weight of the shadow band 
at large k is rather small in Fig. 9. The finite spectral 
weight between the holon and spinon excitation peaks is 
due to the finite system size.^^ 

Taking V/t = 6, we can study the spectral function 
across the Hubbard, non-Hubbard and dominant 4/cf 
regimes with increasing ^. The results are shown in 
Fig. 10, and reveal that the signatures of spin-charge 
separation are fully preserved. Whereas the holon disper- 
sion reflects the noticeable increase of the charge velocity 
with increasing ^, see Fig. 7, the spinon excitations re- 
main virtually unchanged by the interaction range, again 
in accordance with the results for Saiq^oj) in Fig. 7. The 
spectral weight of the shadow band is significantly en- 
hanced compared to V/t = 1 (Fig. 9). On approaching 
the strong-coupling region at larger ^, the upper Hub- 
bard band (visible in the insets of Fig. 10) becomes al- 
most completely flat. Similar to Fig. 10, a gap is visible 
at /cf in Fig. 10(c) [and also in (b) but much smaller]; we 
have verified that this gap is a finite-size effect. 

Our findings in the metallic region of the phase di- 
agram are consistent with the experimentally observed 
coexistence of a small Kp (implying extended-range in- 
teractions) with signatures of spin-charge separation in 
photoemission measurements; a good example is TTF- 
TCNQ.^^ On the other hand, the finite interaction range 
does not provide an explanation for the possible ab- 
sence of clear spin-charge separation in self-organized 
gold chains. ^'^^ We comment on the latter case in the 
conclusions. 

Finally, we show in Fig. 11 the single-particle spectrum 
in the insulating CDW phase at V/t = 9 and ^ = 10. 
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FIG. 10. (Color online) Single-particle spectral function 
A{k,(jj — fi) for n = 0.5, F/t = 6 and different screening 
lengths ^ from exact diagonalization with L = 20. Insets: full 
energy range, revealing the upper Hubbard band. 



The dynamical charge and spin structure factors for these 
parameters were presented in Fig. 8. We find a charge 
gap (equal to 0.2(1) in the thermodynamic limit), and 
backfolded shadow bands related to the 4/cf charge order 
which are visible in the inset of Fig. 11. The spectrum 
appears to evolve continuously across the CDW tran- 
sition. In particular, the holon band is well visible in 
Fig. 11, whereas it has been found to separate into two 
domain walls for much stronger Coulomb interaction.^^ 
The single-particle spectrum of a quarter-filled CDW 
state has also been calculated using the bosonization 
method. In the absence of dimerization, no singularities 
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FIG. 11. (Color online) As in Fig. 9 but for V/t = 9 and ^ = 
10, corresponding to the insulating CDW phase, see Fig. 2. 
The inset shows a logarithmic density plot of the spectrum, 
revealing backfolded shadow bands related to the 4/cf charge 
order. 



exist near (note that in our numerical calculations, we 
cannot distinguish between singularities and excitation 
peaks of finite width). The spectrum may also depend 
on the details of the interaction potential. 



V. CONCLUSIONS 

In this work, we have studied the effects of the electron- 
electron interaction range in one dimension using exact 
numerical methods. We have obtained the Luttinger 
liquid interaction parameter Kp as a function of the 
Coulomb matrix element V and the screening length ^ 
which, in combination with Luttinger liquid theory, de- 
fines the phase diagram of the model. In addition to the 
Hubbard regime 1 > Kp > 1/2, we have explored the 
non-Hubbard regime Kp < 1/2, the case Kp < 1/3 with 
dominant 4/cf charge correlations, and the insulating 
CDW state which exists at quarter filling for Kp < 1/4. 
We identified an important length scale l/Zcp for Kp] Kp 
strongly depends on the screening length for ^ < 1//cf, 
whereas it decays very slowly for ^ l/kp. Our results 
indicate that the lattice model with a finite (but pos- 
sibly large) interaction range can be described by Lut- 
tinger liquid theory if higher-order umklapp terms are 
taken into account. This case is therefore distinct from 
the unscreened 1/r potential which falls outside the Lut- 
tinger liquid description. ^^'^^ For the unscreened poten- 
tial, numerical results suggest the existence of a metal- 
lic quasi Wigner crystal state with Kp = 0.^^ For our 
choice of a screened Coulomb potential, which is both 
convex and monotonically decreasing with increasing dis- 
tance, Kp always decreases with increasing interaction 
strength or range, as compared to enhanced metallic be- 
havior observed in extended Hubbard models as a result 
of competing nearest-neighbor and next-nearest neighbor 
interactions. Interestingly, the small values of Kp ^1/4 



10 



observed in recent experiments on gold chains, as well 
as previously in quantum wires, carbon nanotubes and 
quasi- ID materials, can only be achieved for large values 
of the interaction strength and/or range. 

We have calculated the static and dynamical charge 
and spin correlation functions, and found good agreement 
with the expectations based on Luttinger liquid theory. 
Upon decreasing Kp by increasing V and/or ^, charge 
correlations become strongly enhanced, reminiscent of al- 
though not identical to the quasi Wigner crystal. Our 
results for the real-space density- density correlations are 
consistent with Luttinger liquid behavior on length scales 
beyond the screening length and deviations on smaller 
length scales. 

The 4A:f correlations lead to a pronounced Bragg 
peak in the dynamical density structure factor. The in- 
teraction range strongly modifies the velocity of long- 
wavelength charge excitations, whereas the spin velocity 
only depends on the onsite repulsion. Throughout the 
Luttinger liquid phase, spin-charge separation is clearly 
visible in the single-particle spectrum. Finally, in the 
insulating charge- density- wave phase, we observe back- 
folded shadow bands. 

An important question to be addressed in future 
work, motivated by experiments on self-organized gold 
chains, ^'^^ is the impact of spin incoherence on the spinon 
and holon signatures in photoemission spectra. The en- 
ergy scales for low-energy charge and spin excitations are 
determined by the corresponding velocities Vp and V(j. As 
explicitly shown in this work, Vp increases with increas- 
ing ^, whereas Va- does not depend on the interaction 
range. Therefore, the charge and spin energy scales can 
be well separated for sufficiently large ^. In the regime 
'^p ^ the 2/cf spin correlations can be suppressed at 
finite temperatures, whereas the charge sector remains 
coherent. This scenario may explain the rather inco- 
herent angle-resolved spectrum of gold chains, which at 
the same time show clean LL power-law behavior in the 
angle-integrated density of states.^ 
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APPENDIX: CTQMC 

The general formulation of the weak-coupling CTQMC 
method allows the simulation of problems with long- 



range interactions in imaginary time and/or space. ^^'^^'^^ 
Retarded interactions (i.e., nonlocal in time), which es- 
sentially correspond to the electron-phonon problem, 
have been considered in Refs. 52-54. In this appendix, 
we provide technical details for the application of the 
CTQMC method to a Hamiltonian of the form (1). 

Although such simulations are in principle straight- 
forward, we have encountered difficulties which are ul- 
timately related to the strong-coupling character of the 
problem considered in this paper. The algorithm is quite 
similar to the case of electron-phonon interactions, and 
has been implemented both at finite temperatures and at 
T = (with a projection parameter 0).^^ 

Our starting point is Eq. (1), which we write as 



H = ^ e{k)hk + ^ P{r) {hi - n) {hi^r - n) . 
k 



(A.IO) 

Here n is the average density, the interaction accounts 
for fluctuations around the paramagnetic saddle point, 
and P{r) is a probability distribution; we also defined 
V = '^rV{r). During the simulation, vertices corre- 
sponding to interactions over a distance r are proposed 
with probability P(r) = V{r)/V. 

To circumvent the negative sign problem, and following 
Ref. 34, we rewrite the interaction as 



^ P(r) (f 



'i-\-ra' 



(A.ll) 

Here we have introduced an Ising variable 5 = ±1. Up 
to a constant, Eq. (A.ll) is equivalent to the original 
interaction. To avoid the sign problem for y > we 
have the condition n/2 -\- S > 1. 

The average expansion order, which determines the 
computer time, can be evaluated within the finite tem- 
perature approach, giving 



(M) = PVL 



^P(r)((no 



n) {hr — n)) 



(A.12) 

The fact that the form (A.ll) is beneficial for the simula- 
tions at quarter filling and for rather strong interactions 
conffims an empirically derived rule. In order to obtain 
optimal results away from half filling, it is often useful to 
increase the value of S at the expense of a larger average 
expansion order. With the above formulation, we were 
able to extend the parameter regime of applicability for 
the weak-coupling CTQMC method, and exemplary re- 
sults are shown in Fig. 6. However, the strong-coupling 
regime remains out of reach. 
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